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We theoretically investigate the interplay between charge ordering and magnetic states in quasi- 
one-dimensional molecular conductors TMTTF2X, motivated by the observation of a complex vari- 
ation of competing and/or coexisting phases. We show that the ferroelectric- type charge order in- 
creases two-dimensional antiferromagnetic spin correlation, whereas in the one-dimensional regime 
two different spin-Peierls states are stabilized. By using first-principles band calculations for the 
estimation for the transfer integrals and comparing our results with the experiments, we identify 
the controlling parameters in the experimental phase diagram to be not only the interchain transfer 
integrals but also the amplitude of the charge order. 
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Low-dimensional molecular conductors provide a fruit- 
ful stage to study strong electron correlations in the pres- 
ence of large quantum fluctuations and coupling to lattice 
degrees of freedom [1]. The observed phase transitions in- 
volving spin, charge, and lattice degrees of freedom are 
summarized in the form of pressure (P)-temperature (T) 
phase diagrams for different families. It is common to 
stabilize a different ground state even for relatively small 
pressure variations. Presumably, these changes are trig- 
gered by small variations in lattice constants for a given 
material while maintaining the same geometry of con- 
stituent molecules at room T (isostructural) [2]. A gen- 
eral goal is to identify the parameters controlling the 
ground states and trends in the nature of the elementary 
excitations. An example of interest is the TMTTF2X 
(TMTTF = tetramethyltetrathiafulvalene; X, monova- 
lent anion) family of molecular solids [Sj, in which the 
tuning of charge order, by way of applying pressure, ap- 
pears to play a role in controlling the magnetic states [4| . 
Specifically, decreasing the charge order amplitude by the 
use of applied pressure is associated with the antiferro- 
magnetic (AFM) transition Tn ^ 0, and clarifying the 
relevant physics for this behavior is of interest in the field 
of quantum magnetism. 

In the quasi-one-dimensional family of TM2X (TM: 
TMTTF or tetramethyl-tetraselenafulvalene=TMTSF), 
the key parameter has been widely accepted as the di- 
mensionality {D) tuned by the relative increase of inter- 
chain transfer integrals by P [l, [HI, Figure [T] shows 
the recently updated phase diagram [l, 0, 0, [s| • Amaz- 
ingly, a wide variety of phase transitions appear by ap- 
plying P or a replacement of X (chemical P). The phases 
latest revealed are in the left side where a ferroelectric- 
type charge ordering (FCO) transition was found 0, Isl^ 
it has a strongly correlated nature [lo|, leading to mag- 
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FIG. 1: (Color Online) A schematic phase diagram for TM 
salts [3. The ambient pressures for TMTTF2SbF6 and 
TMTTF2PF6 are shown. DM, FCO, AFM, SP, and SC 
represent dimer-Mott, ferroelectric charge ordering, antifer- 
romagnetic, spin-Peierls, and superconducting states, respec- 
tively. The dashed line is a crossover while the solid lines are 
phase transitions. The inset shows the arrangement of TM 
molecules in the conducting plane. 



netic transitions at low T. Prior to the discovery of the 
FCO phases, the P-D correspondence was invoked to de- 
scribe the phases and transitions appearing at higher P: 
Transport experiments indicate that the correlation gap 
is reduced with P, interpreted as driven by the transverse 
hopping process jhjj , and the system shows a dimen- 
sional crossover [S]. The low-T spin-Peierls (SP) state 
(SPi in Fig. [1]) is destabilized in favor of an AFM state 
(AFMi) [12|, which is consistent with the increase in the 
transverse spin-exchange couplings [13|, \\A \ . 

Difficulties come about when one attempts to apply 
the relation to the left side, the region with FCO. First, 
a discrepancy is easily seen since another AFM phase 
(AFM2) appears at the lowest P and turns into the SP 
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phase (SP2) by applying P, opposite to the SPi ^ AFMi 
variation; the AFM2 and SP2 states both coexist with 
FCO. Another point is that the FCO and AFM2 transi- 
tions behave cooperatively; namely, their transition tem- 
peratures both develop at low P 0, ITU . This is peculiar 
in the sense that, in general, charge ordering tends to de- 
crease the effective spin-exchange couplings lol, isl and, 
therefore, would diminish transitions subject to spin cou- 
plings; in fact, the SP2 phase shows such behavior 0, Isl? 
which is reproduced in theoretical works 

In this Letter, we theoretically elucidate the origin for 
such variations which apparently do not fit to the conven- 
tional practice. Starting by evaluation of transfer inte- 
grals using first-principles band calculations (FPBC), we 
then investigate the roles of electronic correlation and 
electron-lattice couplings on the basis of the effective 
quarter-filled extended Hubbard model (EHM). We will 
show that the complex sequence of phases observed ex- 
perimentally can be reproduced naturally when we add 
the intersite Coulomb interaction as another essential pa- 
rameter, in addition to the interchain transfer integrals. 

The inclusion of the newly found phases in the low-P 
side of the phase diagram was proposed based on NMR 
measurements 0, [H, and the continuous connection of 
phases has been shown by different experiments [lol - 
22|. First we address this in terms of electronic struc- 
ture. We calculate the electronic band dispersions for 
two (TMTTF)2X members situated in the FCO region, 
X=PF6 (SP2 phase) and X=SbF6 (AFM2 phase), within 
FPBC by using the computational code QMAS (Quan- 
tum MAterials Simulator) 23[ based on the projector 
augmented- wave method [2^ with the generalized gradi- 
ent approximation (iij . 



By tight binding fitting to the electron 
bands near the Fermi level we obtain the val- 
ues of transfer integrals in the unit of meV as 
{tai.ta2,h, t,2}={-155, -203, 26.2, -1.31, -3.29} 
for the former and {-149,-207,16.4,-16.4,-9.73} 
for the latter salt (notations are shown in the inset in 
Fig. [1]). The absolute values of tai and ta2 are about 
10 times larger than the other transfer integrals: Both 
salts form a quasi-one-dimensional electronic structure 
along the a axis with dimerization. A measure for the 
D effect is given by \ti)/ta2\ whose values are given as 
0.129 for the PFe salt and 0.080 for the SbFe salt. From 
this point, as far as the transfer integrals are concerned, 
the 2-D in the PFe salt is indeed higher than that in 
the SbFe salt. This is consistent with the semiempirical 
extended Hiickel calculations as well as considerations 
based on their crystal structures (26| . 

Next we investigate the role of Coulomb repulsions on 
top of such an electronic structure, by considering the 
quasi-one-dimensional EHM at quarter-filling in terms of 



holes. The Hamiltonian is given by 



(1) 



where tij is the transfer integral between the neighbor- 
ing sites denoted by (ij), cj^ (Qct) is the creation (an- 
nihilation) operator of a hole on the ith site with spin 
a =t or I, and rii = rii^ + with n^^ = cl^Cia. U 
and Vij are the on-site and the intersite Coulomb inter- 
actions, respectively. From the results of FPBC, we here- 
after set the transfer integrals as tai = —0.8, ta2 = —1, 
and tqi = tq2 = [27] and choose the interchain trans- 
fer integral as a parameter, as inferred from the re- 
sults above. We choose the on-site Coulomb interaction 
to be a typical value for this class of materials [10|, as 
= 4 (~ 1 eV ) and impose a constraint on Vij as 
Vai = Va2 = Vqi = Vq2 = V and Vb = to realize the 
FCO pattern observed in experiments (see Fig. [3j). 

Numerical exact diagonalization on a 4x4 
sites cluster under periodic boundary conditions 
is performed, where we introduce interdimer 
or intradimer charge and spin structure factors 
given by C±{q) = N'^ Z^^^{nfnf)e'^<-^--^^ and 
S±{q) = N-'^,^.{mfmf)e'^<-^--^\ respectively, 
where is the total number of dimers and de- 
notes the center position of the ith dimer. Here, the 
interdimer (+) [intradimer( — )] correlations are detected 
by the summation [difference] in charge and spin den- 
sities within each dimer. 



{n2i 



± ^2i+i)/2 and 

= (m2i ±m2i+i)/2 with rrii = rii^ — respectively, 
where the even (odd) number is labeled as the site for 
the left (right) side in a dimer. 

Figure [2] shows CUq) and S±{q) for tb = 0.1, at F = 
(a) and V = 2 (b)[28|. At F = 0, there is no pronounced 
peak in C±{q); the system is in the dimer-Mott (DM) 
insulating state, where the intrachain dimerization to- 
gether with U leads to a Mott insulator [l^. The en- 
hanced 5'+(7r,g5) and featureless behavior in S-{q) in- 
dicate that the AFM correlation is developed between 
dimers, but only in the a direction due to the 1-D. On 
the other hand, at V = 2, C-{q) has a clear peak at 
q = (0,0), i.e., the FCO correlation, and S±{q) both 
have peaks at gr = (7r,7r). This shows that the develop- 
ment of FCO due to the intersite Coulomb interaction 
induces the two-dimensional AFM correlation between 
the charge rich sites. Note that this happens in spite of 
the fact that the transfer integrals are unchanged from 
Fig. El (a). 

The emergence of 2-D AFM correlation is balanced by 
the degree of FCO and the interchain transfer integral, 
as seen from Fig. [2] (c), where C- (0, 0) and 5'+(7r, tt) are 
plotted on the (h^V) plane. C_(0,0) sharply develops 




FIG. 2: (Color Online) Interdimer and intradimer charge and 
spin structure factors C±{q) and S±{q) at (a) V = and (b) 
y = 2, for [/ = 4 and h = 0.1. (c) C-(0,0) (background 
color) and 5'+(7r,7r) (contour) on the {tb,V) plane: The two- 
dimensional AFM correlation is developed by both tb and V. 
(d) The leading spin-exchange process from the fourth-order 
perturbation in the presence of FCO. 
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FIG. 3: (Color Online) (a) V dependence of the order param- 
eters in the charge densities and lattice distortions for U — 4, 
tb = 0.1, Ki = 0.8, and K2 = I. (b) Ground state phase 
diagram on {tb,V) plane for [7 = 4, Ki = 0.8, and K2 — 1. 
"2DAFM" stands for the enhanced two-dimensional AFM cor- 
relation while the other abbreviations (see text) represents or- 
dered phases. The proposed trajectory for the pressure axis 
in the phase diagram Fig. [1] is shown by the arrow. 



with increasing V while 5'+(7r,7r) increases with increas- 
ing ^5, expected from the interchain spin exchange. The 
noticeable point is that the S'+(7r,7r) peak is rapidly de- 
veloped at large V when FCO is stabilized: The FCO 
state assists in stabilizing the AFM state. As a result, 
the one-dimensional regime in the spin sector is limited 
to the region where both parameters ^5 and V are small. 

The origin of the magnetic properties seen in Fig. [2] (c) 
can be understood by a simplified strong-coupling anal- 
ysis estimating the leading terms of the spin-exchange 
coupling by perturbation calculations with respect to the 
transfer integrals. In the DM state, the spin-exchange 
coupling between dimers along the a axis is simply given 
by J a = —tai/Ud^ while that along the b axis is given 
by J5 = —Atl/Ud-, where Ud is the effective on-site 
Coulomb interaction for the dimer units [lo|- Thus, the 
2-D is enlarged toward \tai\ ~ 2^5, namely, ^5 ~ 0.4; 
in fact 6'+(7r,7r) shows a maximum around this value 



for V = (the large % region is not shown). On the 
other hand, in the basis of the FCO state in the limit of 
large U and V, the charge localizes on every other site 
along the intrachain a axis and on the nearest-neighbor 
sites along the b axis (see Fig. [2] (d)); then the spin- 
exchange coupling between these sites for the former is 
given by Ja ~ —^'^1,i^1l2I{^^^'^) fourth-order 
perturbation whose spin-exchange process is shown in 

Fig. [2] (d) [29], while for the latter J5 r ^l/U from 

the second-order perturbation. Although ^5 is small com- 
pared to tai and ta2, Jh can become the same order com- 
pared to Ja due to the effect of V . Then the 2-D in the 
magnetic state increases and the AFM state is induced. 

It is known that electron-lattice couplings invoke the 
SP transition and various types of phase transitions with 



lattice modulations within one-dimensional models 



17, 



Here we investigate such effects in our quasi- 
one-dimensional system by considering both Peierls- and 
Holstein-type electron-lattice interactions coupled to the 
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EHM, given in addition to Eq. ^ as 

Hp = - ^ Uju,j{clcj^ + H.c.) + ^ E (2) 

{ij)a,cr {ij)a 

where Uij and are the renormahzed lattice distortions 
treated here as classical values and the corresponding 
spring constants are given by Ki and i^2, respectively. 
The Peierls distortions are treated along the a axis [writ- 
ten as {ij)a pairs in Eq. (j2j)]: namely, only intrachain cou- 
plings are considered, to account for the one-dimensional 
quantum effects. Using the Hellman-Feynman theorem 
under the constraint '^ij = 0, we can obtain Uij 

and Vi (= (rii) / K2) self-consistently by using the ground 
state expectation values for bond operators and charge 
densities [13, IH, Hoj . In the following, we show results for 
an 8 X 2 site cluster [si^ at iCi = 0.8 and K2 = I, under 
the antiperiodic and periodic boundary conditions along 
the a and b axis, respectively. From the four-lattice peri- 
odicity modulations along the a axis, we can define order 
parameters by following Ref. [l^ as, for the FCO state, 
n^^: for the coexistence of FCO and SP tetrameriza- 
tion (FCO+SP), {'^co+sp5 ^co+sp}- state 
without FCO (DM+SP), {^i:3M+sP,^DM+sp} [see Fig. 
[3] (b) for schematic representations] . Because of intrin- 
sic dimerization (tai7^a2), two fold lattice distortion u^^^ 
always exists. 

Figure [3] (a) shows the results for ^5 = 0.1 (the same as 
Figs. [2] (a) and (b)), as a function of V. As V increases, 
first a phase transition occurs as DM-hSP FCO+SP 
due to the effect of V [l^ and then to the FCO state with- 
out SP tetramerization; the SP state becomes unstable 
by the development of the two-dimensional AFM corre- 
lation controlled by the FCO that we have seen above. 
The ground state phase diagram on the (^5, V) plane 
is shown in Fig. [3] (b). The two kinds of SP states are 
suppressed with increasing ^5, due to the increase in the 
interchain spin exchange, while, as we have seen in Fig. [3] 
(a), V also diminishes the SP states. We confirm that 
both S±{q) have sharp peaks at gr = (±7r, ±7r) in the 
FCO state, while only S-^{q) has peaks at = (±7r, ±7r) 
in the DM state: These tendencies are the same for the 
case without electron-phonon couplings shown in Fig. [2l 

Based on the above results, we finally discuss our re- 
sults in relation to the complex variation of phases in 
the phase diagram in Fig. [TJ To establish a correspon- 
dence between ^5 in our calculation in Fig. [3] and our 
estimations based on FPBC mentioned above, we need 
to divide the latter by 2 due to the 8x2 cluster hav- 
ing a "ladder" geometry; then for the PFe and SbFe 
salts, this gives ^5 = 0.065 and ^5 = 0.040. By con- 
sidering the experimental ground states for the PFe salt 
(FCO-hSP) and the SbFe salt (FC0+2DAFM), we can 
deduce that they are positioned as indicated in Fig. [3] 



(b), where the SbFe salt has larger value of V. This 
is consistent with the fact that the SbFe salt has larger 
transfer integrals along the diagonal ql and q2 bonds, 
namely, larger overlap between the molecular orbit als, 
which results in larger values of intersite Coulomb re- 
pulsions [111, Vqi and Vq2^ favoring the FCO pattern. 
The smooth evolution of phases with applied P suggests 
that the system follows along the arrow in Fig. [3] (b). 
Specifically, with applied P, transfer integrals reflect- 
ing the overlap between the molecular orbitals are more 
sensitive compared to the inter-site Coulomb repulsions, 
which are approximately a function of intermolecular dis- 
tance [33 • Then, the variation of ground state with P is 
now given by FC0+2DAF (AFM2) ^ FCO+SP (SP2) 
^ DM+SP(SPi) ^ DM+2DAF (AFMi) states, which 
agrees with the variation in Fig. [TJ As for the case of 
chemical P, namely, with the variation among different 
X other than PFe and SbFe, our work suggests that 
a careful reconsideration for each compound should be 
made for how to allocate "ambient P" positions, where 
the anisotropic parameters sensitively reflect the ground 
state; we leave them as a future problem. Our results in- 
dicate that the dimensional crossover in magnetic states 
is controlled by not only but also V; inducing the FCO 
state is essential to understand the sequence of phase 
transitions in TMTTF salts. The apparently confusing 
cooperative behavior in the FCO and AFM2 states is 
naturally understood based on our scenario. 
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